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It is shown theoretically how to directly obtain the energy distribution of a weakly interacting gas 
of bosons confined in an optical lattice in the tight-binding limit. This is accomplished by adding 
a linear potential to a suitably prepared lattice, and allowing the gas to evolve under the influence 
of the total potential. After a prescribed time, a spectral transform is effected where each (highly 
non-local) energy state is transformed into a distinct site of the lattice, thus allowing the energy 
distribution to be (non-destructively) imaged in real space. Evolving for twice the time returns the 
atoms to their initial state. The results suggest efficient methods to both measure the temperature 
in situ, as well as to cool atoms within the lattice: after applying the spectral transform one simply 
needs to remove atoms from all but a few lattice sites. Using exact numerical calculations, the 
effects of interactions and errors in the application of the lattice are examined. 
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I. INTRODUCTION 

Despite the strong interest in the properties of ultra- 
cold bosons confined in optical lattices [3, H, H, 0, [E] , the 
accurate determination and efficient control of the tem- 
perature remains a persistent problem. The initial load- 
ing of the atoms into optical lattices can cause significant 
heating if it is not performed adiabatically Q which be- 
comes increasing difficult as the lattice deepens 0] ■ Like- 
wise, adiabatic changes to the lattice potential itself can 
both increase or decrease the temperature d, [j| [l(| • In- 
strinsic problems such as laser motion [TlT |. atomic mo- 
tion [12] , or impurities [HI, can also adversely affect 
the temperature. The inability to achieve very low tem- 
peratures in optical lattices is one of the primary ob- 
stacles to achieving and controlling interesting quantum 
states in these systems [ill EE E3 • 

There is currently no efficient way to directly measure 
the temperature of the bosons when quantum mechanical 
tunneling within the lattice is appreciable. This weakly 
correlated regime is applicable for sufficiently shallow lat- 
tices, weak interactions, or large numbers of atoms at 
each site [l8|. The standard method is to release the 
atoms from the trap, and is fully destructive (l9| . The 
shadow-image of the fully expanded cloud then corre- 
sponds approximately to the momentum distribution in 
the trap, though the presence of atomic interactions mod- 
ifies this simple picture somewhat [2(| [2l[. In contrast, 
in the weak tunneling (strongly correlated) limit there 
are several proposed methods for inferring the temper- 
ature, including measuring the density distribution of 
doubly occupied sites (22| and the atomic number fluc- 
tuations [23l | . 

Unfortunately, directly cooling weakly correlated 
atoms in an optical lattice is not straightforward. Only 
one procedure has been proposed, based on a combina- 
tion of sophisticated quantum optics and many body ef- 
fects [24|. In contrast, under conditions of weak tun- 
neling several techniques have been proposed [25| . The 



standard evaporative cooling technique for bosons in har- 
monic traps [26|] is predicated on the correspondence be- 
tween the energy and spatial extent of the single-particle 
quantum states. Atoms far from the trap center tend to 
have a higher energy, so their gradual removal combined 
with rethermalization can rapidly cool the particles. This 
approach unfortunately cannot be applied to optical lat- 
tices even if they are shallow, however, because all single- 
particle energy states span the length of the lattice. 

This work proposes that optical lattices can be suit- 
ably prepared in order to allow for the direct measure- 
ment of the temperature of a weakly correlated Bose gas, 
and also to efficiently cool it. The main idea is that if 
one could vary the depth of the optical lattice quadrati- 
cally along its length, then to an excellent approximation 
the low-energy single-particle states of the system are the 
discrete-space analogs of the harmonic oscillator Hermite 
polynomials, known as Krawtchouk functions [27],[28j]. As 
discussed in detail below, this could be accomplished us- 
ing at least two distinct and readily available methods. 
First, one could change the focal length of one of the op- 
tical lattice lasers, so that the intensity maxima of the 
two beams are displaced relative to one another. Second, 
one could add an additional optical lattice at a small an- 
gle relative to the first. In both cases, if the temperature 
is already relatively cold, then an in-situ non-destructive 
image will reveal a bimodal density distribution much 
like those found for bosons in harmonic traps [2j| . Evap- 
orative cooling can then be performed directly 

For larger tem pera tures where the tight-binding band 
is close to full [20] and all single-particle states are 
strongly overlapping, the energy distribution of the gas 
can be obtained by performing a full spectral transfor- 
mation. This is accomplished by adding a small linear 
potential (3lj . and allowing the system to evolve in time. 
Single-particle states with low energies will evolve with 
high probability into the occupation of lattice sites on 
the left, while high-energy states will migrate to sites on 
the right. Thus, after a prescribed evolution time the 
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spatial profile of the atomic density yields the initial en- 
ergy distribution. Further cooling can then be effected 
simply by preferentially removing atoms toward the right 
by imaging with a resonant laser. Time-evolving further 
performs the reverse transformation, mapping the single- 
particle states back to their original representation. 

While this approach is exact in the limit of non- 
interacting particles (and therefore much in the spirit of 
evaporative cooling ideas in harmonic traps) and only 
for a particular choice of tight-binding Hamiltonian, it is 
robust with respect to both atomic interactions and to 
errors in the application of the optical lattice, as shown 
in detail toward the end of this manuscript. It is equally 
applicable to one, two, or three dimensional systems. 

The manuscript is organized as follows. The basic 
formalism for atoms in optical lattices and the tight- 
binding approximation are reviewed in Sec.|TTJ The con- 
ditions under which it is possible in principle to perform 
a spectral transform are discussed in Section IIIII The 
Krawtchouk functions are introduced in Section IIV1 and 
the condition under which these may effect a spectral 
transform is also shown. The robustness of the results to 
the presence of interactions and approximations to the 
optical lattice potential are explored in Section [Vj and 
the results are discussed and summarized in Section [VTl 



II. ATOMS IN OPTICAL LATTICES 

Optical lattices are standing waves formed by overlap- 
ping, focused lasers. In the case of a one-dimensional 
(ID) lattice, the lasers can be considered to be counter- 
propagating. Assuming that the lattice is oriented along 
the axial z direction and is axisymmetric, the lattice po- 
tential is [32| 



V\(r, z) = -V(z) exp(-r 2 /ro) sin 2 (fcz), 



(1) 



where r is the transverse coordinate, V(z) > is the 
lattice depth that may depend on the axial coordinate, 
ro is the transverse beam radius, and k = 2n/X is the 
wavenumber of the light with wavelength A. This ex- 
pression assumes that the confinement is strongest (the 
potential is most negative) at points of largest laser in- 
tensity (i.e. when the frequency c/A of the laser is smaller 
than that of the main atomic resonance, known as red- 
detuning [26[). Allowing for the possibility of an addi- 
tional external potential V e (z) aligned along the lattice 
axis, the Hamiltonian for interacting bosons is 



H 



d 3 rip^ (r, z) 



-—\7 2 + V l {r,z) + V c (z) 
Am 



tp(r, z) 



d 3 rd 3 r'^(r, z)ip j (r' , z')ij}(r', z')^(r, z), (2) 



where ip^(r,z) creates a boson field at position (r, z), 
g = 4:Trh 2 a/m is the two-body coupling constant [26j . 
and a is the s-wave scattering length. For Rb, for ex- 
ample, the value of the scattering length is a w lOOap = 



5.29 nm [33j, with ao the Bohr radius. It is convenient to 
express all energies in terms of the atomic recoil energy 
Eu = h 2 k 2 /2m = 2tt 2 H 2 /m\ 2 . The coupling constant is 
then g = (2\ 2 a/ir)E R . 

In the limit of very strong transverse confinement ^34] , 
all of the boson fields will occupy the lowest-energy eigen- 
state of the confining potential in the transverse direc- 
tion. Neglecting the possible z-dependence of the lat- 
tice potential for the time being, one obtains V\(r) « 
— V(l — r 2 /r 2 ,) = — V + Vr 2 jr 2 which is approximately 
harmonic with effective frequency u> — \/2V jmr^ and 
harmonic oscillator length I = ^h/mtu. Making the sub- 
stitution i/>t(r, z) fn y / 2/n£ 2 exp(-r 2 /2£ 2 )^ (z) and inte- 
grating over the radial coordinate gives a rescaled cou- 
pling constant gin — g/n£ 2 — Anahw = 4:W(aX/r a )E R , 
where V — V /Er is the lattice depth in recoil energies. 

Focusing the optical lattice lasers generally leads to 
an axially-dependent lattice depth V(z) that varies over 
the Rayleigh length. As a result, the lattice depth is 
usually largest at the center and exponentially decreasing 
on either side V{z) = Vq exp(— z 2 /zq), though with a very 
long beam waist zn- Beam waists for optical potentials 
as long as 125 /im [35[ and for optical lattices as small as 
60 /im (36j have been employed in recent experiments. 

An excellent description of ultracold atoms in optical 
lattice potentials of the moderate depth discussed above 
is provided by the Bose-Hubbard Hamiltonian [37] . The 
Bose field operators are expanded ip(z) = biw(z ~ Zi) 
in a basis of Wannier functions w(z — z{) that arc highly 
localized near a lattice site i. Assuming that the ener- 
gies related to the particle dynamics and the temperature 
UbT are both small compared with the bandgap, one ob- 
tains the Bose-Hubbard Hamiltonian 381: 



H 



(ij) 



U 



^K-lH^&fo., (3) 



where (ij) denotes a sum over nearest neighbours (i.e. 
j = i ± 1), b\ creates a boson at site i, and n, = b\b i 
is the density operator for site i. For sufficiently deep 
lattices Vb > 5-Er, one need only consider tunneling of 
atoms between nearest-neighbor sites and interactions 
of atoms within the same lattice site; the amplitudes 
for hopping between next-nearest-neighbor sites and the 
strength of interactions between atoms on neighboring 
sites are smaller by at least an order of magnitude, re- 
spectively [37l | . 

If the external potential V e (z) is slowly varying over 
the lattice spacing, the hopping coefficient is 



Jij = / dzw*(z — Zi) 



h 2 9 



w(z-zj). (4) 



In the very deep lattice limit where the value of the lat- 
tice depth at site i (in recoil energies) is large Vi = 
Vi/En 3> 1, the Wannier functions can be approximated 
by the ground-state Gaussian wavefunctions for each sep- 
arate well w(z) w 



l 2d l^fd^ with d 
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Then the hopping amplitudes dU) are given by Jij /Er = 
{■K 2 /A)Vi exp(— tt 2 \f~V~i / 4) . A better approximation for 
the Wannier functions instead gives [f| [22j 



Ji 



Er 



16—3/4 

— V l exp 

7T 



-2\jV 

Meanwhile, the on-site interaction term is 
U = ,9id / dz\w(z)\ 4 . 



(5) 



(6) 



Again approximating the Wannier functions by harmonic 
oscillator functions, one obtains 



U 



a\k— 3/4 



/32^F 3/4 . 



7r r r 
Finally, the site-dependent local energy is 



(7) 



dzw*(z — Zi) 



-^V 2 + V l (z) + V c (z) 
2m 



w(z- 



(8) 

This latter term is a constant energy offset for homoge- 
neous lattice potentials and no exernal potential, and is 
usually dropped. In the present case, however, the lattice 
amplitude will vary slowly over the lattice spacing A/2. 
This will give rise to a spatial variation of Si, which for 
large lattices is equivalent to a smoothly varying external 
potential even in the absence of any V e (z). Again using 
the local harmonic oscillator solutions for the Wannier 
functions, one obtains Ei/Er ss \fv~i- 

The onset of strong correlations for bosons in quasi- 
1D optical lattices with an average of (n) ^S> 1 atoms per 
site occurs when U/J > 2.2(n) [a, yjj, |4(|. Combining 
Eqs. ([5|) and |7|), one obtains the criterion 



ro 



> 2.2(n) 



(9) 



The value of ro is set by imposing the quasi- ID criterion 
huj > VEr, or r /A <C (irW)^ 1 . This condition is in 
fact much more restrictive than necessary because at low 
temperatures T cx J (for example, the BEC transition 
temperature for a 3D lattice is fcsT c = 5.59J [5]) only 
the first band is occupied, and the bandwidth in the ID 
tight binding limit is w = 2 J < V. With Hlu = VE R , 
a = 5.29 nm, and A = 800 nm, and a typical value (n) = 
1000 one obtains V^max < 16 and ro « 32 nm to ensure 
that the atoms remain in the weakly correlated regime. 
This transverse length scale corresponds to a confining 
frequency lu/2tt w 100 kHz, a typical value in current 
experiments With (n) = 10 000 instead, one obtains 

^ max ~ 25. 

Far in the weakly correlated regime, the amplitude for 
nearest-neighbor tunneling from site to site tends to be 
larger (and therefore more important) than the strength 
of particle interactions. As a first approximation, we may 
therefore set U = 0; the consequences of U ^ will be 



explored in Sec. [V] Furthermore, it will be assumed that 
V e (z) can be chosen in such a way to ensure that is 
approximately constant, so that it can be ignored (this 
assumption will be further justified in Sec. IV A[) . The 
resulting non-interacting tight-binding Hamiltonian can 
then be written 



H — J i,3 
{ij) 



i)(3\, 



(10) 



where the site state kets are defined by \i) = b\\0) and 
|0) is the vacuum state. 

In principle, the position indices i and j are unre- 
stricted, but in practice the lattice has a finite number 
of occupied sites. This could arise naturally because of 
the inherent spatial variation of the lattice potential dis- 
cussed above, in which larger number of atoms will oc- 
cupy the central region where the potential is deepest. 
To enforce occupation of a specific part of the lattice, 
one can apply 'end-cap' potentials formed from tightly 
focused blue-detuned lasers [13, [3f|. In this case, one 
can effectively enforce the 'hard-wall' boundary condi- 
tions i e {0, N} with j = i ± 1 when < i < N, but 
j = 1 only for i = and j = N — 1 only when i = N . 



III. SPECTRAL TRANSFORM 



Introduction 



As is well-known from elementary quantum mechan- 
ics, an arbitrary quantum state \ip(j)} defined on N + 1 
discrete sites labelled by index j can always be expressed 
in terms of an arbitrary basis of dimension N + 1 



N 



(11) 



71 = 



where the c n are generally complex coefficients. The set 
of states {\f n {j)),n £ 0,N} is assumed to be both or- 
thonormal 



N 



^('PmU^fnij)) = (VmlVn) = <W (12) 
j=0 



and complete 



(13) 



where the Kronecker delta is defined as 6j k — 1 only if 
j = k and is zero otherwise. Multiplying both sides of 
Eq. pip on the left by (f m (j)\ ar "l summing over j, one 
obtains the coefficients using (| 1 2p : 



(14) 
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One could have obtained this directly simply by inserting 
the complete set, \ip) = J2 n \<Pn) (<Pn\4>} ■ 

The relevance of the above discussion to the present 
work is that the coefficients c„ can be considered as the 
representation of the original wavefunction in the 'spec- 
tral basis' rather than the original site basis. While the 
elements of the original wavefunction in the site basis are 

(M) = -00') = ^2c n (j\fn) = ^2c n tp n (j), (15) 

n n 

the elements of the same wavefunction in the spectral 
basis are 

3 3 

Eqs (fT5|) and (p~6|) are transforms of each other through 
the basis functions <p n (j)', the second can be immedi- 
ately obtained from the first from their orthonormal- 
ity (fi"2")) , and the converse follows from their complete- 
ness (fl"3)) . These generalize the Fourier transforms that 
are specific to the regular lattice case, in which case 
tp n (j) oc exp[i2mij / (N + 1)]. 

The purpose of the current work is to show that the 
coefficients c„ in Eq. (fTT|) can be obtained directly in 
experiments without having to explicitly carry out the 
transformation (|16p . In particular, if one could imple- 
ment the spectral transform 

l^-e^ln) (17) 

with arbitrary phases <p n , then all particle densities would 
be transformed as |0(,7')| 2 — > \cj\ 2 : 

\m\ 2 = mm = J2 c ^\^Y, c *n(fm\j) 

n m 

- ^CnC^e^-^mWj) = \ Cj \\ (18) 

m,n 

Thus after the spectral transformation, the magnitude of 
the coefficients in the expansion (jTTj) would be obtained 
simply from the values of the density at each lattice point. 

B. Transform by Hamiltonian Evolution 

Consider a single particle hopping on N+l sites of a ID 
lattice with hard- wall boundary conditions, as defined by 
the non-interacting tight-binding Hamiltonian pTJ|) with 
Jy = Jjj and J—i,o — Jjv,iv+i = 0. This Hamiltonian 
has N + l eigenvectors \vj) with eigenvalues ej. Now 
consider the perturbed Hamiltonian H = Hq + Hi with 
completely local (i.e. on-site) potential 

N 

#i=£^|j>01 (19) 

3=0 

with real Vi . (Note that the Vi is an external local poten- 
tial that is applied in addition to the V e (z) chosen above 



to ensure constant e^). The objective is to determine 
the Vj for a given choice of such that time-evolution 
governed by the full Hamiltonian H effects the desired 
spectral transformation (|17p . but now with \v n ) taking 
the place of the \<j> n ) above: 

exp(-iHt a /h)\vj)=e^\j), (20) 

where t s is some fixed time satisfying the condition for 
any j S {0, N} and <pj is a possible arbitrary phase. 

It is insightful to write the desired unitary U = 
exp(—iHt s /h) in its spectral decomposition 

N 

u = Y,^M-^' k ts/h)\v' k )(v' k \, (21) 

fc=0 

where e' k and \v' k ) are the eigenvalues and eigenvectors of 
H, respectively. The transform (|20[) can then be written 

N 

U\ Vj ) =Y,eM-iekts/h)\v' k )(v' k \v j ) = e»>\j). (22) 
Clearly, if 

(v' k \v j )=eMKts/K)e i ^(v' k \j) (23) 
then one would guarantee that 

U\v s )=e^ \j) = e**'\j) (24) 

\k=0 I 

where the term in parentheses is unity because the eigen- 
basis of H is complete. 

It would be preferable if the system would evolve back 
to the original eigenbasis after an additional time t s , at 
least up to an overall phase. This is equivalent to requir- 
ing that 

N 

U 2 = ^M-^ k t s /h)\v' k ){v' k \ = Ie ie , (25) 

k=0 

where / is the identity matrix of dimension N + l and 9 
is some phase. This condition is most generally obtained 
when 2e' k t s /h — 2im — 6 for any k, n. Suppose that we 
choose n — k m with m some arbitrary integer (m = 1 
implies the spectrum is linear, m = 2 gives a quadratic 
spectrum, etc.). Then e' k = e' + ck m , with e' Q = —h8/2t s 
and c = %-K/t s . Eq. (|22"|) can now be written 

N 

U\vj) = e^-'^^^WhWM (26) 

k=0 

N 

= e*/ a 5>-** fcm K>Kh) 

fc=0 
N 

= e i9 / 2 ^(-l) fe |^)KK)=e^|i). 

fc=0 
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Thus, in order to perform a spectral transform that is 
invertible by repetition, one requires 

(v' k \ Vj ) = (-i)V<*'- fl / 3 ><«;|j) Vk,je{0,N}. (27) 

In short, the overlap of the eigenvectors of the full Hamil- 
tonian H with those of the unperturbed Hamiltonian Hq 
must be trivially related to their overlap with the site 
basis vectors. Furthermore, the perturbation to Hq must 
be local (i.e. be strictly an externally applied spatially 
dependent potential) , and the spectrum of the full Hamil- 
tonian must be an integer power in the quantum number. 

Consider the homogeneous optical lattice with = ,/ 
on N sites. Hard-wall boundary conditions force all 
functions to zero at the sites labeled by i = —1 and 
i = N + 1 . The unperturbed eigenvalues in ascending or- 
der are e k = -2 Jcos[7r(fc + 1)/(N + 2)] with k e {0,iV}, 
and the eigenvectors are similar to Chebyshev polyno- 
mials, \v k ) = y/2/{N + 2) sin [n(k + + l)/(N + 2)], 
where k is the eigenvalue index and j £ {0, N} corre- 
sponds to the lattice site. It is straightforward to verify 
that choosing Vk = £fe, i.e. where the perturbation cor- 
responds to the eigenvalues of the original unperturbed 
Hamiltonian in ascending order, indeed yields perturbed 
eigenvectors \v' k ) that satisfy the desired overlap relation 
(v' k \vj) = (— l) k+ i {v' k \j) . Unfortunately, the spectrum of 
the perturbed Hamiltonian is not of the desired special 
form e' k = e' Q + ck m . Conversely, if a perturbation is cho- 
sen to enforce the eigenvalue criterion, then the overlap 
criterion cannot be satisfied. The unfortunate conclusion 
is that no perturbation to a uniform optical lattice can 
be chosen in order to effect a spectral transform. The 
remainder of this work is devoted to showing that both 
of these conditions can be nevertheless be simultaneously 
satisfied with atoms in specially prepared optical lattices. 



IV. KRAWTCHOUK LATTICE 



Consider the site-dependent tunneling amplitudes 



J, 



Jy/(i + l)(N-i), 



(28) 



where i £ {0, N — 1} and J is a bare tunneling coeffi- 
cient with units of energy. The associated weighted lat- 
tice Hamiltonian has been recently inve stig ated in the 
context of quantum walks [U, |42|, HH, Hj, the quan- 
tum mechanical generalization of the classical random 
walk [i^, |46[. In particular, with these coefficients in the 
governing Hamiltonian, a quantum state initialized in the 
left-most lattice site will evolve with perfect probability 
to the rightmost lattice site after a fixed amount of time. 
For this reason, this Hamiltonian and its generalizations 
have been denoted a perfect quantum wires [471 ] . Prop- 
agating further for the same amount of time returns the 
state to the first site. Likewise, states at any site will be 
transferred to sites symmetrically located around the lat- 
tice center, so time-evolution effects a perfect reflection 
of any state about the midpoint. Much of the interest in 



this and related systems stems from their potential uses 
for quantum computation and communication, because 
of their ability to perfectly distribute quantum states and 
entanglement [48l l49l] . 

In the present context the Hamiltonian (jTUJ) with tun- 
neling coefficients ([28]) is useful because it possesses both 
a linear spectrum and exhibits periodic behavior. Indeed, 
the eigenfunctions turn out to be Krawtchouk functions, 
which are the discrete-space analogs of the Hermite poly- 
nomials. Furthermore, it will be shown below that per- 
turbing the quantum wire by a linear potential will yield 
a full Hamiltonian whose states satisfy the criterion (|2T[) . 
thereby yielding a full spectral transformation between 
energy and coordinate space. 



A. Unperturbed Hamiltonian 

The spectrum of the unperturbed Hamiltonian (|10p 
is obtained by solving Ho\v k ) = efc|ufc) Expanding the 
eigenfunctions in the site basis \v k ) = afj \j) with 

(k) . rr> 

dj some unknown expansion coefficients, one can write 
the eigenvalue equation as 

a f ( J <+i,<|i + vw) + j *-mI* - 1> w> 

-e k \j))=0, (29) 

which can be conveniently rewritten as the three-term 
recurrence relation 



t (fe) , t ( fe ) O) n 

J j-j-l a ]-l + J j,j+l a j+l - e kdj - U. 

In terms of the quantum wire tunneling amplitudes 
this becomes 



(30) 



-e k a\ = 0, 



(31) 



where the eigenvalues €j = tj/ J are now dimensionless. 
The recurrence relation (1311) defines the normalized 



Krawtchouk polynomials ip^^z) [27l. [2 



^ Np + ^- 2p)j - k ^\k) = o, 



vWi-p) 



(32) 



where the Krawtchouk polynomials can be expressed in 
terms of either hypergeometric or Jacobi functions: 



I ^UF,(-j,-*;-«;ij, 



6 



1-p 



(-l)'(l-p) 



(N-j-k,-N-l) 



n/2 im-jV- 




k\(N-k)\ 
(33) 



xP, 



with p G R a free parameter. Comparison of Eqs. (f3T[) 
and (|32p immediately requires p = 1/2 to eliminate the 
j (position) dependence of the last coefficient. The spec- 
trum of the unperturbed Hamiltonian is therefore linear, 
ej. = — N + 2k. The eigenvectors are 



Mi)> 



2 W / 2 V fc!(iV-fc)! J 



,(JV-j-fc,-JV-l) 



(3), (34) 



where the position-dependence is now included explicitly. 

In the limit of large N and small k, the Krawtchouk 
polynomials reduce to Hermite polynomials [U H(| , re- 
flecting the fact that the Hamiltonian for the quantum 
wire is the discretized version of the harmonic oscilla- 
tor. This is easy to verify directly from the recurrence 
relation (|3Tj) in the limit N — > oo and j <C N: 

,(fc) 



fjNaf^ + VU + W4% +( N ~ 2k ) a f ] = °> ( 35 ) 

which is equivalent to the recurrence relation for the nor- 
malized Hermite polynomials Hj(x) p5l| 



y/j&j-l + y/j + lH j+ i ~ xV2Hj = 0. 



(36) 



Evidently, x — \/2k neglecting the overall constant term. 
The effective trapping (oscillation) frequency is therefore 
U! = 2J/h. 

B. Full Hamiltonian 

The j-dependence of the third coefficient appearing in 
the recurrence relation (f3"2"| indicates that adding a linear 
potential Vj oc j in (|19[) to the unperturbed Hamiltonian 
would still yield eigenstates that are Krawtchouk func- 
tions. In this case the potential would need to be of the 
form 



v - 1 ~ 2p ■ 

yielding the softened linear spectrum 
1 



(-Np + k) 



(37) 



(38) 



and the full expression (|33| for the eigenvectors \v' k (j)) = 

It remains to be shown that there is a choice of p for the 
full Hamiltonian that can satisfy the main condition (f27|) . 
In fact, as will be proven below, there is indeed a choice 
of p that ensures 



(v t \v„ 



-1)* 



(39) 



which effects the spectral transformation U\v m ) — 
(— l) m \m) for all m, and another choice that performs 
the mirror-symmetric version U\v m ) = {—l) m \N — m). 
In terms of the Krawtchouk functions, condition (|39[) can 
be written as 



E 

3=0 V J 



2^/2 I m J i 1 ' 



P 



1-p 



J/2 



N \ p {N-j-k,-N-l) ( 2 _ 
k)n [p 



P \ ^ ( N\ p(N-m-k-N-l) 

i-p) m \p 

(40) 

Eq. P0|) can be thought of as an (N + l)-dimensional 
matrix equation K\Ki = K3, with the sum over j 
on the left-hand side effecting the product of the two 
Krawtchouk matrices [52|, [53l 



K 



(AO 



(~l) fe 
2 N/2 



N 



P 



(N-m-j,-N-l) 



(3), (41) 



K. 



(N) 



= ("If 



1-p 



J/2 



x p(N-j-k-N-l) ( 2 __ l 
3 \P 



(42) 



which are clearly symmetric under the interchange of row 
and column indices j <-» m or j <-> k respectively, and 
the right hand side given by the symmetric Krawtchouk 
matrix 



K 



(AT) 



V 



1-p 



1/2 



p{N~ 



-k,-N-l) 



The factor 



in the sum (l40l) serves as 



- 1 

(43) 
'met- 



ric tensor,' as described in greater detail below. The 
Krawtchouk matrix elements are found by fixing N and 
treating the two variables in the superscript of the Jacobi 
polynomials as the row and column indices, respectively. 
Thus, the columns of the symmetric Krawtchouk matri- 
ces K[ N ^ and correspond to the eigenvalue-ordered 
eigenvectors of the unperturbed and perturbed Hamilto- 
nians Ho and H\ on N + 1 sites, respectively. 

As has been recently discussed in Ref. HH, Krawtchouk 
matrices correspond to the 'Hamming contraction' of 
tensor products of identical matrices. The hierarchy 
of Krawtchouk matrices is mathematically expressed as 

ifW = r |(if( 1 ))® Ar |, where is the Nth order 

Krawtchouk matrix of dimension iV + 1, and r is an oper- 
ator effecting the Hamming contraction. The contraction 
is effected by labeling the rows / columns of the contracted 
matrix by the Hamming weight of the row/column index 
(written in binary) for the original matrix, and summing 
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all matrix elements whose row (column) indices share 
the same Hamming weight. For example, the first-order 

Krawtchouk matrix is the Hadamard operator H 
' I 1 1 



K 



1 

71 



i 



i 



(44) 



(2) 

and the second-order Krawtchouk matrix iVJ is the con- 
traction of the tensor product of two Hadamard opera- 
tors: 



K 



(2) 





(45) 



1 



Equation (|4H| can therefore be interpreted as the 
product of uncontracted tensor products of first-order 
Krawtchouk matrices, 



K. 



(1) 



which can be alternatively expressed as 



(46) 



(47) 



The 'weight factor' 



appearing in the sum of 



Eq. (|4"0")) can now clearly be understood as accounting 
for the multiplicities for each Hamming weight on the 
integers between and 2^, and therefore serves as a 
'metric' for matrix-matrix multiplications in contracted 
space. Thus, to prove the equality (|4U|) for arbitrary N 
one need only find the value of p satisfying the N = 1 
case: 




which yields the unique solution 



V2 



(48) 



(49) 



It is easy to verify that this solution indeed satisfies 
Eq. (|40|) for arbitrary choices of N. Inserting this result 
into Eqs. ([57)) and gives the linear potential 



V 3 = 2 j 



(50) 



and the resulting linear spectrum for the perturbed 
Hamiltonian 



v / 2(-iV + 2fc). 



(51) 



Because c = 2y2 in Eq. (|26|) . the time to effect the spec- 
tral transform is 



7T fh 



(52) 



Eqs. (f2"8"| . (|5"0|) and (|52"|) constitute the main results of 
this work. 

To perform the mirror-symmetric spectral transform 
U\v m ) = | TV — to) one needs to obtain the value of p that 
ensures 



(v'k\ v m) = (-1 
In this case, one obtains 



, k+m 



v' k \N-m), 



(53) 



K. 



(N) 



= (-If 



k-\-m 



P \ f N 



l-p 



(m-fc,_JV_l) (2 



(54) 



Following the same procedure as discussed above, it is 
straightforward to show that to satisfy condition (I53|) one 
requires p = j (2 + y/2) . This yields the same prefactor 
as the solution above, namely l/y/p(l — p) — 2\[2. The 
spectrum is therefore the same, but the applied potential 
now has the opposite sign, Vj = —2j. Thus, to perform 
an inverted spectral transform where the ground state is 
transformed to the rightmost site rather than the left- 
most, one simply requires the same linear potential but 
with a negative slope. 



C. Discussion of Main Results 

One might be surprised that neither the applied linear 
potential (|50p nor the required time (|52|) depend explic- 
itly on the number of lattice sites. Clearly, as N grows, 
the initial wavepacket must traverse an increasing num- 
ber of lattice sites during the implementation of the spec- 
tral transform. It is important to recall, however, that 
the actual tunneling amplitudes (|2"5)) do indeed increase 
with N. That said, it is not obvious that their values are 
just what are required to ensure that t s remains constant. 
Consider the time it takes on average to tunnel between 
adjacent lattice sites, U^+i w The total time 

required for a wavepacket to propagate from one end of 
the lattice to the other is then approximately 



T : 



J 



N 1 

V- 1 — 



(55) 



In the large- TV limit, one obtains the asymptotic result 
limjv^oo T = (H/J)tt = h/2J. This is exactly the period 
for one oscillation of the original Krawtchouk oscillator: 
the eigenvalues = 2k are equivalent to the frequency 
2J/h = 2tt/T, which yields T = h/2J. 

The fact that the Krawtchouk matrices represent con- 
tractions of tensor products of identical matrices has 
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in fact been anticipated in the literature, though in a 
completely different context. The quantum wire tun- 
neling amplitudes ([28]) are obtained by 'collapsing' the 
hypercube graph to a weighted linear graph, by render- 
ing indistinguishable vertex labels with the same Ham- 
ming weight [4lj |. The hypercube is constructed by ten- 
sor (Cartesian) products of two-site linear graphs. The 
eigenvectors are therefore tensor products of the two so- 
lutions for the two-site graph. These solutions are rep- 
resented by the Hadamard matrix H. The Hamiltonian 
and its eigenvectors for the quantum wires (the collapsed 
hypercubes) are therefore contractions of the associated 
hypercube quantities. 

It should be pointed out that the relationship between 
Krawtchouk functions and spectral transforms has been 
discussed previously in the literature, primarily in the 
context of optical and digital filtering [54], [55|. These 
ideas were extended to quantum mechanical systems, 
where the Krawtchouk transform can be considered as a 
fractional Fourier transform • That said, the idea pre- 
sented in this paper of performing the spectral tranform 
by time-evolution under the influence of a local operator 
is entirely new. 

It is also important to mention that the discrete nature 
of the Hamiltonian is a crucial ingredient that enables the 
spectral transform to be performed simply by turning 
on a linear potential. Suppose we attempted the same 
approach in the continuum limit N — > oo in which the 
Krawtchouk functions reduce to Hermite polynomials. In 
this limit, the inhomogeneous tight- binding Hamiltonian 
is replaced by a harmonic oscillator Hamiltonian with a 
potential V(z) = with uj — 2J/h as discussed at 

the end of Sec. IIV Al Adding a linear potential V(z) = 
2J[z/ (X/2)] = 4Jz/X, the continuum limit of the discrete 
linear potential Vj — 2j required to effect the spectral 
transform, the total potential becomes 



V(z) 



■raw 



z + 



4J 
rriLU 2 \ 



8J 2 



; 2 A 2 ' 



(56) 



The linear potential simply shifts the origin of the effec- 
tive oscillator potential to the left, so the eigenstates of 
the unperturbed Hamiltonian are instantaneously placed 
at the classical turning points of the shifted potential. 
The subsequent evolution is simple harmonic motion, 
with all wavepackets initially propagating to the left. Un- 
der no circumstances will certain eigenstates propagate 
to the right! 

One more note should be made before discussing the 
feasibility of performing this spectral transform in ac- 
tual optical lattices in the next section. The potential 
that must to be applied (V(j) — 2Jj) in order to effect 
the spectral transform corresponds exactly to the eigen- 
values of the original (unperturbed) Hamiltonian, listed 
along the diagonal in increasing order. Curiously, this 
was exactly the same criterion found to satisfy the con- 
dition (f2"T)) for the uniform lattice, discussed at the end of 
Sec. IIII1 though in that case the full spectral transform 



could not be implemented because the resulting eigen- 
value spectrum was not an integer power-law. This ob- 
servation could provide a tantalizing clue to generalizing 
the spectral transform by Hamiltonian evolution for other 
discrete systems. 



V. APPROXIMATIONS AND INTERACTIONS 

The implementation of a spectral transform by Hamil- 
tonian evolution described in the previous section as- 
sumed that there is only one particle in the entire lat- 
tice. This is of course not generally the case in real ex- 
periments, where there are many atoms and collisions 
between pairs of particles is not negligible. It is therefore 
important to investigate the efficacy of the approach in- 
cluding the effects of particle interactions. This is accom- 
plished by including the interaction term proportional to 
U in the Bose- Hubbard Hamiltonian ([3]) , which has been 
so far neglected. 

Likewise, it has been assumed that the site-dependent 
tunneling coefficients Ji,i+i in (|28p could be implemented 
in actual experiments. In practice, the tunneling ampli- 
tudes can be made to vary continuously over the length 
of the lattice, as discussed in Sec. [Hi but this variation is 
unlikely to yield exactly the correct values of the Jj,i+i. 
This section considers separately the robustness of the 
spectral transform when the exact conditions discussed 
in Sec. Mil are no longer satisfied. 



A. Quadratic Optical Lattice 



The exact tunneling coefficients (|2"5)l required for the 
spectral transform would not be straightforward to im- 
plement experimentally. As discussed in Sec. [HI a natural 
variation of the lattice depth arises from the laser beam 
waists zq, V(z) = Vb exp(— z 2 /z 2 ). With the substitu- 
tion z — (A/2)i, the axial variation for the discrete case 
is Vi = Vq exp[— (A/2zo) 2 i 2 ]. Inserting this into the tight- 
binding expression for the hopping coefficients J5]) and 
assuming that zq 3> A one obtains 



M+i 



Jo 



1 



A 

2z~ 



where Jn = 



16 F 3/ V2 



(57) 



and terms of order 



(A/2zo) 4 have been neglected. Assuming zq = 100 jxm, 
A = 0.8 //m, and Vq = 30, this quadratic approxima- 
tion gives Ji,i+i/En fts J (l + 7.5 x 10^ 5 i 2 ). Can such 
a quadratic variation of the tunneling coefficients repro- 
duce the desired values ([28]) ? 

In the vicinity of the lattice center i — (N — l)/2 + j 
with j a small integer j <C N, the tunneling coefficients 
are indeed approximately quadratic: 



= J 




(58) 
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N + l 



J 1- 



2f 



(A + I) : 



(59) 



neglecting terms of order 0(j 4 ). While this is an excel- 
lent approximation near the lattice center, the tunneling 
coefficients decrease more quickly near the lattice edges; 
clearly, if j = (N — l)/2 in the expression above one ob- 
tains J(n~i)/2,(n+i)/2 ~ N/4, which is much larger than 
the exact tunneling coefficient ^(jv-1)/2,(jv+i)/2 = VN- 
An alternative approximation that reproduces the de- 
sired values at both j = and j = ±(N — l)/2 would 
be 



N + l 



J 1 



4 f 



(a + i)(va + i) 2 



(60) 



Comparison of Eqs. (|57|) and (|60|) reveals a problem, 
however: the coefficients of the site-dependent term have 
opposite signs. This is because the actual tunneling co- 
efficients for a focused laser beam ((57)) are smallest when 
the lattice depth (laser intensity) is greatest. In contrast, 
the desired coefficients (|60|) are supposed to be largest at 
the center, so that the lattice depth must be shallowest 
at this point. 



1. Axially displacing lattice beams 

In order to prepare an optical lattice with the desired 
properties, one could axially displace the lattice beams 
relative to the each other, so that the intensity maxima 
are located a distance d ^> zo apart. The lattice poten- 
tial would then be V{z) = ^V {e- Z ^ z o + e -(^-d) 2 /^). 
Written in terms of the new variable x = z — d/2 centered 
midway between the two maxima, the potential becomes 



V(x) = V d/2 e 



2/2 

~ x ,z ° cosh 



xd 



(61) 



where V 



d/2 



Voe~ d / Az ° is the zeroth-order potential 
depth at the lattice center equidistant from the centers 
of the two beam maxima. Expanding about x — gives 
an approximate quadratic envelope 



V{x) « V m + V, 



d/2 



2z 2 



(62) 



corresponding to tunneling coefficients 



J(x) « J, 



d/2 



1 - 


x 2 


fd 2 




'24 


U 


-i) 



Vd/2- 



where J d/2 = J^-V^J^e" 2 ^ v d ' 2 . Comparison of 



(63) 



and (|6"3"j) immediately gives J = 2J d / 2 /(N + 1) and 



A 2 



V d/2 - 



32 



(N + 1)(VN + 1) 2 ' 

(64) 



where the substitution x = (A/2)j has been made. Of 
course, to enforce the hard-wall boundary condition, end 
caps [3, Hll must also be applied transversally at the 
optical lattice amplitude maxima located at z — and 
z = d. 

According to Eq. ([SJ, the spatial variation of the 
potential depth (|62p would yield a site-dependent en- 
ergy shift that is not included in the assumed hop- 
ping Hamiltonian (1 10[) . Using the harmonic oscil- 
lator approximation to the Wannier functions yields 

e{x)/Er w JV d / 2 x 2 d 2 /4zq ignoring the constant off- 



set energy yV d / 2 . It is therefore necessary to manually 
apply an inverted parabolic potential, centered at the 
depth minimum, that can serve to cancel this effective 
quadratic potential energy contribution. 

Let's consider what these results imply for actual 
experiments, assuming Rb atoms confined in a red- 
detuned optical lattice with wavelength A = 0.8 fim. 
With Vq = 50Er and zq — 100 /im, one could obtain 
a lattice with N — 201 sites choosing d w 255 /xm. 
The potential depth at the midpoint between beam max- 
ima is V d / 2 w 10 En, deeply within the tight-binding 
limit. The coupling constant is then calculated to be 
J ps 2.3 x 10~ i Efi, though the actual coupling constants 
vary from 0.01-Er at the edge to 0.02Eft at the midpoint 
between laser maxima. An external potential must be 

applied at all times, V e (a;) « — 5.1(x/zq) 2 so that the 
magnitude of the external potential reaches almost one 
recoil energy at the lattice edges. Applying the addi- 
tional linear potential V^(x) — 4J(x/X) then allows for 
the implementation of the spectral transform in a time 
t„ = hn/2V2J fa 210 ms. 

Likewise, for a lattice with N = 301 sites and Vq — 
20Er, one requires d ~ 223 /im which yields V d / 2 ~ 
5.8Er. The external quadratic potential that needs to be 

applied is V° e (x) » -3(x/z ) 2 - With J = 4.6 x fO" 4 ^ 
(the actual coupling constants vary between 0.036-Er and 
0.069Er), the time required is faster than for the exam- 
ple above, t s = 107 ms. It might at first glance appear 
paradoxical that the N = 301 case should be faster than 
the N = 201 case, but recall that the total time t s is set 
by the parameter J, which for the latter case is twice as 
large as the value for the former case. 

In order to verify that the quadratic optical lattice is 
indeed able to perform the spectral transform, an ex- 
plicit numerical simulation of the time evolution was per- 
formed. The governing tight-binding Hamiltonian is the 
unperturbed Hq defined in Eq. (|I0|) , but with the approx- 
imate tunneling coefficients Jij = J(j) = J(x) given in 
Eq. (|6"3"|) ; its eigenvalues and eigenvalues and eigenvectors 
are and \vk), respectively. For simplicity, the initial 
state was chosen to be a Boltzmann distribution 



P 



= y^^cxp(-e k /k B T)\v k )(v k \, 



(65) 



where p is the thermal density matrix, T is the temper- 
ature in Kelvin, and ks ~ 1-38 x I0 -23 J/K is Boltz- 
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FIG. 1: Numerical results for the spectral transform. The pa- 
rameters for the simulation are ksT — 40 Er, zq = 100 /im, 
Vo = 50 Eh, and N = 100. The left and right graphs corre- 
spond to d = 200 pra and 300 pra, respectively. The results 
for the perfect tunneling coefficients (|28[) are given by dashed 
lines, while those based on tunneling coefficients (JSJ) derived 
from the actual lattice potential (|61|) are given by solid lines. 
The insets depict the density profile of the initial thermal 
distribution for the actual lattice. 

mann's constant. The time-evolution of the density ma- 
trix is obtained in principle by numerically integrating 
the equation of motion ihdp/dt = [H + 4Jx/A, p] where 
[, ] denotes the commutator and the correct linear poten- 
tial has been explicitly added. Equivalently, because the 
particles are assumed to have no interactions, one may 
simply allow each eigenvector of the unperturbed Hamil- 
tonian to evolve separately under the influence of the 
Schrodinger equation, and then add their contributions 
to obtain the time-dependent density matrix 

P {t) =J2eM-tk/k B T)\v k (t)}{v k (t)\. (66) 

fc 

The spatially dependent particle density is then read- 
ily obtained using the relation n(j,t) = (j\p(t)\j) or 
n(x,t) = (x\p(t)\x). 

The numerical results for the representative set of pa- 
rameters ksT = 40 Er, z — 100 ^m, Vo — 50 Er, and 
N = 100 are shown in Fig. [TJ for the choices d = 200 pra 
and 300 pm. Two cases were considered. In the first, the 
tunneling coefficients were chosen to be the exact ideal 
tunneling coefficients (|2"5)) . As expected, the initial ther- 
mal distribution over energies is exactly reproduced after 
a time t = it as a spatial distribution over lattice 

sites j, shown as the dashed lines in the figure. In the 
second case, the tunneling coefficients were obtained for 
the actual lattice by inserting the spatially dependent op- 
tical lattice potential (|6ip into the general expression (0 
appropriate in the tight-binding limit. The results are 
shown as solid lines in the figure. It is clear that for 
both values of d, the actual optical lattice does an ex- 



cellent job of reproducing the spectral transform. This 
indicates that a careful adjustment of the separation of 
the intensity maxima for the two optical lattice lasers 
(subject to rf > Zq) is not necessary to ensure a suc- 
cessful implementation of the spectral transform. These 
general conclusions are also found to be independent of 
temperature. 

It is also worthwhile to point out that the ability to 
perform the spectral transform directly in the optical 
lattice greatly simplifies both the measurement of tem- 
perature and the ability to cool the atoms. The initial 
particle density (shown as insets in Fig. [I} is very broad, 
with tails that are strongly truncated by the lattice edges. 
This would make a direct measurement of temperature, 
obtained by fitting the tails to an exponential function, 
quite unreliable. Furthermore, while one could attempt 
to cool the initial cloud directly by ejecting particles far 
from the lattice center, one would have little control over 
the actual energies that were removed from the system. 
In contrast, the particle density after the time propaga- 
tion is almost an exact spatial map of the initial energy 
distribution. The exponential tail (at least for UrT < Vo) 
is well-resolved, allowing an excellent determination of 
temperature. Furthermore, evaporative cooling on this 
state simply requires the removal of all atoms occupy- 
ing sites with index j > 1, followed by another spectral 
transform back to the energy basis. 



2. Spatial variation by interference 

While the approach described above of offsetting the 
focal distances of the two laser beams is able to repro- 
duce the desired spectral transform, many experiments 
are performed under the assumption that the hopping 
amplitudes are approximately constant over the region of 
interest. It would be preferable to employ a homogeneous 
optical lattice (with hard-wall boundary conditions) for 
experiments, and then to only introduce the spatial inho- 
mogeneity when temperature measurements or cooling is 
desired. Unfortunately, changing the focal length in the 
middle of an experiment is not feasible. 

An alternative approach to modulate the optical lat- 
tice is through interference. Suppose that the homo- 
geneous lattice is formed from two counter-propagating 
laser beams each with amplitude A and frequency Q in 
the standard lin-_L-lin configuration, i.e. where the two 
polarizations e perpendicular to one another. The elec- 
tric field is then 

E(z) = Ae x e t27rz/X e~ mt +Ae v e- t27rz/x e- tQt 

= Ae l27rz/X e~ int (e x + e a e-* 47rz/A ) . (67) 

The spatial variation of the polarization vector gives rise 
to a periodic atomic polarizability, which in turn yields 
the lattice potential V(z) cx cos 2 (27rz/A). If one now 
introduces a second beam with electric field amplitude 
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B <C A oriented in the yz-plane at a small angle 9 with 
respect to the z-axis, 



where 



-iSlt J i27r[zcos(0)+i/sin(0)]/A 



+ £ye -* 2Tr [ 2:cos ( e )+asm(e)]/A| 

then the total lattice potential becomes 



V(z) = V Q { cos 



2ttz 



a cos 



2ttz 



cos(0) 



(68) 



. (69) 



where a oc |S| 2 /|^4| 2 <C 1 is the relative depth of the 
second lattice beam. 

For small angle 9 < 2tt so that cos(0) » 1- 2 /2, 
Eq. (jnU) can be rewritten as 



V(z) ps U <^ cos 



2?rz 



a cos 



27TZ 27TZ 



A 



(70) 

where the 'beat wavelength' has been defined as Ab = 
2 A/6* 2 3> A, for reasons that will be clear shortly. The 
potential then becomes 



V(z) 



1 + a cos 



2?rz 



cos 



2?rz 



/ 47TZ \ / 27TZ 

-asmf — )sm( — 



(71) 



If a is small then the last term in the above ex- 
pression is small and can be neglected. The poten- 
tial is therefore described by a spatially varying depth 
Vq [1 + a cos fiirz/Xb)] 2 which is periodic at the beat 
wavelength. Unlike the case a = 1 often encountered 
in physics, however, the potential depth never vanishes: 
the maximum depth V^n ax ps Vb(l + 2a) corresponds to 
zsiO, while the minimum depth for z ps A&/2 = A/0 2 is 

Knin Pa fo(l - 2a). 

Expanding about the minimum z — x + Ab/2 in the 
regime x <C Ab one obtains 



V(a;) « 



1- 2a 



47r 2 a 

\2 ; 

A 6 



COS 



(72) 



where as usual V = V/Er. Thus, again the lattice depth 
is approximately a quadratic function of the lattice po- 
sition. This leads to a quadratic site-dependent energy 
contribution JSJ: 



e(x) 



' Fp 27r 2 ax 2 
1 - 2a A 2 



(73) 



neglecting site-independent energies. Substituting this 
expression into Eq. ([5]) yields the site-dependent tunnel- 
ing coefficients in vicinity of the depth minimum 



J{x) « J 



Af,/2 



4a7r^ ( y/Vo - - 



(74) 



>^A b /2 — ^0 



1 + 2a 



(75) 



is the tunneling amplitude at the lattice minimum in 
terms of the bare tunneling amplitude Jo- 

Comparison with Eq. (foT))) gives J = 2J Xb / 2 /(N + 1) 
and 



(7V + l)(v / iV+l) 2 _ "Ub 



U - 



(76) 



It is again useful to insert numbers that would be relevant 
to experiments with ultracold atoms. With 9 — tt/25, 
the beat wavelength is Ab = 126. 6A or approximately 
253 lattice spacings. One should therefore be able to 
effect a spectral transform with N < 201. Choosing Vo 
I (I E R and N = 101, condition ||7BJ yields a ps 0.216, 
giving J ss 9 x 10~ 4 E R and therefore t s « 54 ms. Note 
that an inverted parabolic potential V e (x) ps — 18(ir/Ab) 2 
must be applied at all times to cancel the site-dependent 
energy contribution (|73|) . Likewise, choosing N = 201 
requires a ps 0.058 and t s ps 170 ms. For larger values of 
N, one would require smaller values of 9 to ensure that 
A 6 > N(X/2). 

The main advantage of this approach is that tempera- 
ture measurements and/or cooling can be effected at any 
stage of an experiment with ultracold atoms in optical 
lattices: one simply adds the second lattice beam and the 
appropriate quadratic and linear potentials. Of course, 
these must be added adiabatically in order to avoid heat- 
ing the system in the process. A conservative estimate of 
the adiabaticity timescale would correspond to the time 
taken by a wavepacket to propagate fully across the lat- 
tice a few times, i a diab ~ 2H/J from Eq. (f5"5"]) . For the 
N = 101 example above, one obtains i a diab ~ 100 ms so 
that a ramp-up of the additional beam and the external 
quadratic potential over 200 to 300 ms would ensure that 
the system would be minimally perturbed. 



B. Effect of Interactions 

Until now the interaction term in the Bose-Hubbard 
Hamiltonian ([3]) has been completely neglected. In prin- 
ciple, the strength of particle interactions can be adjusted 
close to zero through the use of Feshbach resonances [56[ , 
but in general the gas will not be perfectly ideal. Un- 
fortunately, a full numerical simulation of the spectral 
transform, as was discussed in the subsection above, is 
not straightforward in the presence of interactions. The 
interactions lead to mixing among the energy levels dur- 
ing the time propagation that are impossible to model in 
an average way without introducing additional assump- 
tions. 

Instead, the spectral transform was numerically simu- 
lated by exact propagation of the many particle Hamil- 
tonian, but for a smaller number of lattice sites N = 10. 
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FIG. 2: Fidelity F as a function of the number of particles 
N p obtained by exact numerical simulation of the spectral 
transform. With the initial state the ground state of the un- 
perturbed interacting Hamiltonian, the fidelity is defined as 
the density in the first site at time t s = ty/2\/2J. The solid, 
dashed, and dot-dashed lines correspond to U = 0.1J, 0.2J, 
and 0.5J, respectively. 



The tunneling coefficients corresponded to the ideal ex- 
pression ([28]) . and the values of U in Eq. ([3]) and the num- 
ber of particles 1 < N p < 10 were varied independently. 
The largest number N p — 10 for N = 10 corresponds to 
affilbert space dimension of (N+N p -iy./(N p -l)\N\ = 
92 378; the computational time for the simulation of the 
spectral transform became unpleasantly long for larger 
N p . The goal was to determine the fidelity F of the 
spectral transform; i.e. if the initial state was the fcth 
eigenvector of the unperturbed (but interacting) Hamil- 
tonian, then F is the particle density in the fcth site of 
the lattice at the time t s . 

In practice, one need only consider the fidelity of the 
\vq) <-> |0) spectral transformation, because this is rep- 
resentative of all the eigenstates. The results of the nu- 
merical simulation are shown in Fig. [5] When U — 0.1J, 
the fidelity has dropped to only 0.98 when N p = N = 10, 
i.e. at unit lattice filling v = N p /N = 1. The best fit to 
the data indicates that the fidelity drops as a power law, 
F « 1 — 0.023^ 33 . As the on-site interaction strength 
is increased the fidelity drops; when v = 1 one obtains 
F w 0.92 when U = 0.2J and only F = 0.63 when 
U = 0.5J. The best fits in the U = 0.2J and U = 0.5J 
cases are F w 1 - 0.088^ 3 3 and Fal- 0.43z/ 31 , respec- 
tively, both power laws with a similar exponent. Though 
it is impossible to confidently extrapolate to large sys- 
tem sizes with N and N p comparable to those in a future 
experiment, the numerical results suggest that the abil- 
ity to perform a spectral transform to 90% accuracy or 
better would require v < 0.6 and/or U < 0.5J. While 
this value of U /J ensures that the system is in the super- 
fluid regime [5jJ, it is smaller than the value one would 



generally obtain in the ID lattice case (cf. Sec. HI)) . 



VI. CONCLUSIONS 

The scheme described above to perform the spectral 
transform with ultracold atoms has assumed that the 
optical lattice is strictly one-dimensional. While this as- 
sumption has simplified the mathematical analysis and 
the numerical simulations, it is not necessary. In prac- 
tice, many Bravais lattices are generated by two counter- 
propagating laser beams in different (not necessarily per- 
pendicular) spatial directions. One could envisage that 
the lattice beams along each different axis would be 
suitably prepared so that their intensities were varying 
quadratically, using either approach discussed in Sec. [V] 
The tunneling amplitudes would then vary along each 
axis, so that the eigenstates of the lattice (in the non- 
interacting limit) would correspond approximately to the 
product of Krawtchouk functions for each direction. 

Performing a full three-dimensional spectral transform 
would not necessarily require a linear potential applied 
in each direction simultaneously, a task which might be 
difficult in practice. Instead, one could simply apply the 
linear ramps for each dimension in sequence, separately 
cooling the atoms' motions in each spatial degree of free- 
dom. For the cubic lattice for example, eliminating atoms 
from large site values in one direction would not perturb 
the atoms in the perpendicular directions, and cooling 
can be efficiently performed. The presence of interac- 
tions will modify this picture somewhat, but if the inter- 
action strength is kept at sufficiently small (U < J based 
on the one-dimensional simulations) then this should not 
present much of a problem. 

Besides the obvious benefits of temperature measure- 
ment and cooling, the ability to perform a spectral trans- 
formation in situ has other possible future applications. 
For example, the technique could be used to probe the 
excitation spectrum and heating mechanisms of strongly 
interacting systems. At any given moment, the inter- 
actions could be quickly tuned to approximately zero 
through the application of Feshbach resonances, and the 
system's energy spectrum monitored by spectral trans- 
form followed by non-destructive imaging; subsequently 
the system would be returned to its original state by in- 
verting the process. Of course, care would need to be 
taken to ensure that the ramp on and off of the interac- 
tions were slow enough that this process itself would not 
greatly perturb the system. 

As a second example, it is conceivable that the tech- 
nique could be extended to produce a greater variety of 
overall phases than the and it in the present case (|39p . 
In particular, if a suitable choice of tunneling param- 
eters and external potential could be found such that 
the resulting phases would be of the form irkm/N, then 
the spectral transform would be equivalent to a quantum 
Fourier transform, a central operation in several efficient 
quantum algorithms including factoring [58, 59]. In this 
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case, each lattice site (out of a total N — 2") would label 
a possible state in the n-qubit Hilbert space. Though 
inefficient in terms of resources, this notion could pave 
the way toward implementation of interesting quantum 
algorithms using ultracold atoms in optical lattices with 
minimal external controls. 
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